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A reduced Hamiltonian formulation to reproduce the saturated regime of a Single Pass Free 
Electron Laser, around perfect tuning, is here discussed. Asymptotically, N m particles are found 
to organize in a dense cluster, that evolves as an individual massive unit. The remaining particles 
fill the surrounding uniform sea, spanning a finite portion of phase space, approximately delimited 
by the average momenta u+ and U)-. These quantities enter the model as external parameters, 
which can be self-consistently determined within the proposed theoretical framework. To this aim, 
we make use of a statistical mechanics treatment of the Vlasov equation, that governs the initial 
amplification process. Simulations of the reduced dynamics are shown to successfully capture the 
oscillating regime observed within the original iV-body picture. 



I. GENERAL BACKGROUND 



Free-Electron Lasers (FELs) are coherent and tunable radiation sources, which differ from conventional lasers in 
using a relativistic electron beam as their lasing medium, hence the term free-electron. 

The physical mechanism responsible for the light emission and amplification is the interaction between the relativistic 
electron beam, a magnctostatic periodic field generated in the undulator and an optical wave copropagating with the 
electrons. Due to the effect of the magnetic field, the electrons are forced to follow sinusoidal trajectories, thus 
emitting synchrotron radiation. This spontaneous emission is then amplified along the undulator until the laser effect 
is reached. Among different schemes, single-pass high-gain FELs are currently attracting growing interest, as they are 
promising sources of powerful and coherent light in the UV and X ranges. Besides the Self Amplified Spontaneous 
Emission (SASE) setting Q, seeding schemes may be adopted where a small laser signal is injected at the entrance 
of the undulator and guides the subsequent amplification process 0- In the following we shall refer to the latter 
case. Basic features of the system dynamics are successfully captured by a simple one-dimensional Hamiltonian 
model |l5j introduced by Bonifacio and collaborators in 0] . Remarkably, this simplified formulation applies to other 
physical systems, provided a formal translation of the variables involved is performed. As an example, focus on 
kinetic plasma turbulence, e.g. the electron beam-plasma instability. When a weak electron beam is injected into 
a thermal plasma, electrostatic modes at the plasma frequency (Langmuir modes) arc destabilized. The interaction 
of the Langmuir waves and the electrons constituting the beam can be studied in the framework of a self-consistent 
Hamiltonian picture 0], formally equivalent to the one in 0. In a recent paper 0] we established a bridge between 
these two areas of investigation (FEL and plasma), and exploited the connection to derive a reduced Hamiltonian 
model to characterize the saturated dynamics of the laser. According to this scenario, N m particles are trapped in 
the resonance, i.e. experience a bouncing motion in one of the (periodically repeated) potential wells, and form a 
clump that evolves as a single macro-particle localized in space. The remaining particles populate the surrounding 
halo, being almost uniformly distributed in phase space between two sharp boundaries, whose average momentum is 
labeled lo + and uj_. The issue of providing a self-consistent estimate for the external parameters N m , uj + and cj_ is 
addressed and solved in this paper. 

This long-standing problem was first pointed out by Tennyson et al. in the pioneering work 0] and recently revisited 
in 0. A first attempt to calculate N m is made in where a semi-analytical argument is proposed. In this respect, the 
strategy here proposed applies to a large class of phenomena whose dynamics can be modeled within a Hamiltonian 
framework 00] displaying the emergence of collective behaviour 0. 

The paper is organized as follows. In Section [H] we introduce the one-dimensional model of a FEL amplifier 
and review the derivation of the reduced Hamiltonian 0,0]. Section ITTT1 recalls the statistical mechanics approach to 
estimate the saturated laser regime. In Sections II V I to IVD the analytic characterization of N m , 10+ and uj- is given in 
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details and the results are then tested numerically in section lYlII Finally, in Section IVlIII we sum up and draw our 
conclusions. 

II. FROM THE SELF-CONSISTENT iV-BODY HAMILTONIAN TO THE REDUCED FORMULATION 

Under the hypothesis of one-dimensional motion and monochromatic radiation, the steady state dynamics of a 
Single-Pass Free Electron Laser is described by the following set of equations: 

dF = « ' « 

^ = -Ae** - A*e- i9 > , (2) 
dz 
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where z = 2k w pzj 2 / '(7)0 is the rescaled longitudinal coordinate, which plays the role of time. Here, p = 
[a w u> p / (4ck w )] 2 / 3 /j r is the so-called Pierce parameter, (7)0 the mean energy of the electrons at the undulator's 
entrance, k w = 2tt/\ w the wave number of the undulator, uj p = (47re 2 n/m) 1 / 2 the plasma frequency, c the speed 
of light, n the total electron number density, e and m respectively the charge and mass of one electron. Further, 

1 /2 

a w = eB w /(k w mc 2 ), where B w is the rms peak undulator field. Here 7 r = (\ w (l + a 2 u )/2A) is the resonant 
energy, X w and A being respectively the period of the undulator and the wavelength of the radiation field. In- 
troducing the wavenumber k of the FEL radiation, the two canonically conjugated variables are (0,p), defined as 
6 = (k + k w )z — 2Spk w zj 2 /('y}'o an dp = (7— (l)o)/ (p(l)o)- 9 corresponds to the phase of the electrons with respect to 
the ponderomotive wave. The complex amplitude A = A x + iA y represents the scaled field, transversal to z. Finally, 
the detuning parameter is given by S = ((7)0 — 7r)/(2/ 3 7r)i an d measures the average relative deviation from the 
resonance condition. 

The above system of equations (N being the number of electrons) can be derived from the Hamiltonian 



N 2 pf N 

^ = Ey- w + 2 Viv^ sin( ^'" (/ ' ) ' (4) 
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where the intensity / and the phase ip of the wave are given by A = y/l/N exp(— iip). Here the canonically conjugated 
variables are (pj,9j) for 1 < j < N and (J, (p). Besides the "energy" H, the total momentum P = J2j Pj + ^ is also 
conserved. By exploiting these conserved quantities, one can recast the FEL equations of motion in the following 
form for the set of 2N conjugate variables (qj,pj) 



1 N 

= Y,^ qi + 6, (5) 

1=1 




9j 



Pj = ~ 2 VT7 C0S '?^ (6) 



where the dot denotes derivation with respect to z, and qj = 9j — <p mod(27r) is the phase of the j th electron in a 
proper reference frame. The fixed points of system JHJ)-© are determined by imposing qj = pj = and solving: 

1 N 

= Y,sinq t +S = , (7) 



Pj ' 

\/NI 

1=1 



An elliptical fixed point is found for q% = q = 3tt/2. The conjugate momentum solves (p + 5)y / P/N — p + 1 = and 
therefore depends on P/N. We shall return on this issue in the following Sections. 
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For a monokinctic initial beam with velocity resonant with the wave, equations Q , (0) and predict an exponential 
instability and a late oscillating saturation for the amplitude of the radiation field. Numerical simulations fully confirm 
this scenario as displayed in fig. H In the single particle (q,p) space, a dense core of particles is trapped by the wave 
and behaves like a large "macro-particle", that evolves coherently in the resonance. The distances between these 
particles do not grow exponentially fast (as is the case for chaotic motion) but grow at most linearly with time 
(for particles trapped in the resonance with different adiabatic invariants, i.e. essentially different action in the single 
particle pendulum- like description). This linear- in-time departure of the particles appears in the differential rotation 
m fig. [3 while the remaining particles are almost uniformly distributed between two oscillating boundaries. Having 
observed the formation of such structures in the phase-space allowed to derive a simplified Hamiltonian model to 
characterize the asymptotic evolution of the laser 0, [|| . This reduced formulation consists in only four degrees of 
freedom, namely the wave, the macro-particle and the two boundaries delimiting the portion of space occupied by the 
so-called chaotic sea, i.e. the uniform halo surrounding the inner core. 




FIG. 1: Evolution of the radiation intensity as follows from equations Jlj) and ffi) . N = 10 4 electrons are simulated, for 
an initial mono-energetic profile. Here 5 = and 1(0) ~ 0. Particles are initially uniformly distributed in space. 



In [5j we hypothesized the macro-particle to be formed by N m individual massive units, and introduced the variables 
(C, £) to label its position in the phase space. 

The N c = N — N m particles of the surrounding halo arc treated as a continuum with constant phase space 
distribution, f sea (6,p, z) = f c , between two boundaries, namely p+{6, z) and p~(6, z) such that: 

P±=P°±+ P± exp(i8) + p* ± exp(-i6) , (9) 

where p± represents their mean velocity. These assumptions allow to map the original system, after linearizing with 
respect to p± , into Q : 



C = i<S> e K - i<5>* e~ li (10) 

/ ± = -ic 

2 ± 2 

i N. 

2~N N 



\V± = ~u;±V±+i$ (11) 
<i> = i^ V +- V - +i ^e-« + i5Z (12) 



where ^ij 
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A = 
P± = 
P° ± = 



V±/2 



(13) 
(14) 
(15) 



Normalizing the density in the chaotic sea to unity yields f c = l/(27rAw), where Atu := lo + — cj_ represents the 
(average) width of the chaotic sea. The above system can be cast in a Hamiltonian form by introducing new actions 
I± and their conjugate angles ip±: 



V- 



(16) 
(17) 



A pictorial representation of the main quantities involved in the analysis is displayed in fig. [3] In addition: 



$ = iA 



(18) 



The reduced 4-degrees-of-frccdom Hamiltonian reads, up to a constant irrelevant to the evolution equations: 



2N„ 



51 + — — 2a 



/ Ji+ sin(<£ - <p + ) - sfTTZ sin(ip + (p-)] - 2/V7 sm(ip - (), (19) 
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FIG. 2: Phase space portraits for different position along the undulator [z — a) 80, b) 81, c) 83, d) 84/. 27ie differential rotation 
of the macro-particle is clearly displayed. For the parameters choice refer to the caption of Fig. QJ 
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FIG. 3: (q,p) phase space portrait in the deep saturated regime for a monokvnetic initial beam (1(0) ~ 0, Pj(0) = and q 
uniformly distributed in [0, 2ty]). The two solid lines result from a numerical fit performed according to the following strategy. 
First, the particles located close to the outer boundaries are selected and then the expression p±(q) = lo± =p |V±| sin(q + B±) is 
numerically adjusted to interpolate their distribution. Here, \V±\, B± and lo± are free parameters. The numerics are compatible 
with the simplifying assumption B+ = B- ~ 0. 



where 



Nr. 



NAuj 



(20) 



= ^ (21) 



N 

The first four terms represent the kinetic energy of the macro-particle, the oscillation of the wave and the harmonic 
contributions associated to the oscillation of the chaotic sea boundaries. The remaining terms refer to the interaction 
energy. Total momentum is P = £ + I + N c u>. The Hamiltonian (|19|l allows for a simplified description of the late 
nonlinear regime of the instability, provided the three parameters w+, a->_ and N m are given. 

To achieve a complete and satisfying theoretical description we need to provide an argument to self-consistently 
estimate these coefficients. To this end, we shall use the analytical characterization of the asymptotic behavior of the 
laser intensity and beam bunching (a measure of the electrons spatial modulation) obtained in [Tlj with a statistical 
mechanics approach. In the next section these results are shortly reviewed. 



III. STATISTICAL THEORY OF SINGLE-PASS FEL SATURATED REGIME 



As observed in the previous Section, the process of wave amplification occurs in two steps: an initial exponential 
growth followed by a relaxation towards a quasi-stationary state characterized by large oscillations. This regime is 
governed by the Vlasov equati on, rigorously obtained by performing the continuum limit (N — > oo at fixed volume 
and energy per particle) jj, [ill U3 on the discrete system (11l.''>t . Formally, the following Vlasov- wave system is found: 



p-^ + 2(ArCos6> - A y sin 0)-^- , (22) 



% = -P% + 2(A x cose- A y sin 6)?f 
oz oO op 

dA x 

dz 



-5 A v + J fcos6&6dp , (23) 
dA f 

—^-=SA X - fsm8c\9dp . (24) 
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The latter conserves the pseudo-energy per particle 



e= f ^-f(9,p)d9dp-5(A 2 x +A 2 y ) + 2 J (A x sin0 + A y cos9)f(6, p)d6 dp (25) 
and the momentum per particle 



P f(6,p)d6dp+(A 2 X + A 2 y ) . (26) 

A subsequent slow relaxation towards the Boltzmann equilibrium is observed. This is a typical fmite-iV effect 
and occurs on time-scales much longer than the transit trough the undulator 0, 0, 0| . For our calculations we are 
interested in the first saturated state. To estimate analytically the average intensity and bunching parameter in this 
regime we exploit the statistical treatment of the Vlasov equation, presented in In the following, we provide a 
short outline of the strategy. Since the Gibbs ensembles are equivalent for this model, note that the same expressions 
are recovered through a canonical calculation 0, 0, 0] . 

The basic idea is to coarse-grain the microscopic one-particle distribution function f(6,p,z). An entropy is then 
associated to the coarse-grained distribution /, which essentially counts a number of microscopic configurations. 
Neglecting the contribution of the field, since it represents only one degree of freedom within the (N + 1) of the 
Hamiltonian 10} , one assumes 



'(/) = - 



fo fo 



f ode dp 



fo fo 



fodOdp 



(27) 



where the constant fo is related to the initial distribution [l7j . 

The equilibrium is computed by maximizing this entropy, while imposing the dynamical constraints. This corre- 
sponds to solving the constrained variational problem 



S(e,a)= max ( s(f) H(f, A x , A y ) = Ne; P(f, A x , A y ) = Na; / f(9,p)d9dp=l 

f J A x ,Ay 



(28) 



which leads to the equilibrium values 



/ 
.4 



fo: 



-/3(p 2 /2+2A sin0)- \p-n 
e -/3(p 2 /2+2Asin6>)-Ap-^ 



Al 



[35 - A 



S m(9)f(6,p)d6dp, 



(29) 
(30) 



where (3, A and fi are the Lagrange multipliers for the energy, momentum and normalization constraints and, in 
addition, we have assumed the non-restrictive condition ^cos(#i) = Using then the three equations for the 
constraints, the statistical equilibrium calculation is reduced to finding the values of the multipliers as functions of 
energy e and momentum a. These equations lead directly to the estimates of the equilibrium values for the intensity 
/ and bunching parameter b = \ ^exp(iftj)|/iV. 

In the following, we focus on the case of an initially monokinctic beam injected at the wave velocity, while the 
initial wave intensity is negligible, so that e = and a = 0. Moreover we let fo —* oo, which amounts to /x — ► oo in 
cq. (J23). Results are displayed in fig. 0] showing remarkably good agreement between theory and simulations, below 
the critical threshold 8 C — 1.9 that marks the transition between high and low gain regimes. This transition is purely 
dynamical and cannot be reproduced by the statistical calculation. 

Analytically, it turns out that 



fc _ U |3 , , I 1 (2/(3|A|3-2^|)) 
where l n is the modified Bcssel function of order n. In particular for 5 = 0, one finds |j4| 2 = I/N 



(31) 

0.65 and b ~ 0.54. 
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FIG. 4: Comparison between theory (solid and long-dashed lines) and simulations (symbols) for a monoenergetic beam with 
e = 0, a — 0, when varying the detuning S. The dotted vertical line, S = 5 C ~ 1.9, represents the transition from the high-gain 
to low gain regime 0/. 



IV. TOWARDS THE ANALYTICAL CHARACTERIZATION OF N, 



As previously discussed, one can predict the value of the bunching parameter b, using the above statistical mechanics 
description. Clearly, the bunching parameter b depends on the spatial distribution of the particles. From its definition 
it immediately follows: 



b = 



1/2 



sin qi 



(32) 



To proceed we can isolate the contribution relative to the macroparticle from that associated to the chaotic sea. 
We thus obtain: 




Focus on the first two terms of expression l|33|) . We assume the macroparticle to be ideally localized at the 
elliptic fixed point solving Q-©, i-e. set Qj = q = 37r/2, for each individual massive unit belonging to the inner 
agglomeration. Hence, ^2 i^macro cos li = N m cosq — 0. As concerns the second contribution, recalling the expressions 
for the boundaries P-{q) and p+(q) (see caption of fig. |3J), one can formally write: 
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! N c 1 [ 27r f p+{q) , , 

— > cos qi ~ -T7TT1 \ / cosq dpdq 



n c i r 27r 

Nr. 



(cos q)[u+ — W-]dq 



N 2tz{u 

The other contributions in eq. (|33[) can be estimated as follows 



r/ 1 + sing] cos gdg = 0. 

_ UJ -) Jo 



(34) 



- ^ sin% ~-_ (35) 

i^macro 



while 



- ^ sing, ~ — — / (sing) / dpdg 



N ^ N 2tt(uj + -lu) Jr. 

i<£sea v r ' u 

7V r 1 



p-(a) 



TV 2tt(cj + - w_) y 



/ sin g(w+ — w_)dg 
Jo 



V r + | + |7_|)sin^dg = --i— -, (36) 



TV 2-k(lu + — cj-) J Vl +l 1 " 7V2A^' 
where V = \V+\ + \V-\. Inserting into JSBJ yields: 



7V TO 7V C V 

6= — + lV2A^ (3?) 



and finally 



= 2bAuj-V = l _ 2(1 - b)Au 
N 2Auj - V 2Auj - V 

using the relation N m = N — N c . It is worth emphasizing that, neglecting as a first approximation the amplitudes of 
the sinusoidal boundaries, i.e. setting V- = V+ = 0, the previous equation reduces to 

^ = b, (39) 

confirming the relevance of the macroparticle picture to the bunching parameter, a physical quantity of paramount 
importance for the FEL dynamics. Formula i|37|) yields bunching parameter values larger than l|39|l . i.e. implies that 
the sea also contributes to increasing the bunching parameter. The difference increases as the sea gets more populated 
and the width lo + — lu- is reduced (note that our approximations require |V+| + |V_| < lu+ — see fig. |3J). 

V. ESTIMATING THE AVERAGE MOMENTUM OF THE BOUNDARIES w± 

In this section we estimate the unknown quantities u± by characterizing their functional dependence on N m . For 
this purpose we introduce (see schematic layout of fig. |3J) : 

w± = Q ± ^ . (40) 
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FIG. 5: Solid line: Au vs z. Dashed line: lj vs z. Parameters are set as discussed in the caption of fig. QJ 

The problem of estimating u± is obviously equivalent to providing a self-consistent calculation for uj and Aco. The 
latter are both monitored as function of time in fig. [3] and shown to be practically constant. 

In the following, we shall focus on the case of a system which evolved from an initially monokinetic beam and an 
initially infinitesimal wave. It is then convenient to choose the Galilean reference frame moving at the beam initial 
velocity. This translates into the conditions e = and a = 0. The detuning <5 is arbitrary so far. It is of prime interest 
to consider the special case where the beam is injected at the resonant velocity, so that 5 = 0. We shall make this 
additional assumption in section rVIII but the estimates in this section and in the next one do not require it unless 
explicitly stated. 

Consider the conservation of momentum for the original TV-body system Q and focus on the asymptotic dynam- 
ics, which allows one to isolate the contributions respectively associated to the macroparticle and the chaotic sea. 
Averaging over the number of particles yields: 



--TPsea + --J-Pmacro + J = °~0 (41) 

where p sea stands for the average momentum of the chaotic sea and the subscript '0' labels the initial condition. 
To simplify the calculations, we introduced the rescaled intensity J = I/N. As already observed in 0, 0, the 
macroparticle rotates in phase space. This rotation is directly coupled to the oscillations displayed by the laser 
intensity. Averaging over a bounce period z rot , one formally gets: 



^(Psea) + ^(Pmacro) + (J) = (42) 

where (■) stands for the time average. Focus now on (p sea )- Since particles are uniformly filling the chaotic sea, one 
can use the approximation outlined before eq. (see also caption of figEU): 



(Psea) = / [ pf sea (e, P ,z)d6d P )dz=— -// pdpdO- {LJ + +LO_)/2=Q . (43) 



z rot Ji-Q \J J ) 2it{uj + - ui-) 



p-(l) 



As already outlined in the preceding discussion, we assume that the macroparticle oscillates around the fixed point 
and therefore each individual clement constituting the macroparticle verifies the condition (qj) = q. In addition, from 
equation J7J: 
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To proceed in the analysis, we approximate the right hand side in equation 14411 as: 



/ 1 



E sing, 
\Ny/j J- 1 



; sin(cfj) = ■ sing 

V \ / i pmacro V \ / i dm aero 



N m 1 



N VU) 



(45) 



consistently with the argument after (J33J. The above relation is derived by performing a linearization (see Appendix), 
validated numerically and supported a posteriori by the correctness of the results. The contribution of the chaotic 
sea reads: 



777 £ si 



iGsea 



sing. 



-±= E si 



j£sea 



2 71 



p+(g) 



N c 1 

sinfl, ~ 7^=7l — Tk — r / sm <7 / dpdg 



2ir p2ir 

(lu + — uj^ ) sin qdq — / {\V+\ + \V-\) sin 2 gdg 
o Jo 



N c V 
Ny/{J)2Auj' 



(46) 



Merging equations l]44fl. Ij45|) and Ij46(l and recalling (|37(1 and l].'Ufl. one obtains 



{Pmacro) 



N m 1 N c V 1 



V /(J)" TV 2Aw v /p) 



-(J). 



(47) 



Thus the macroparticle moves on the average at the same velocity as the center of the chaotic sea. Inserting equations 
((43(1 and 1)47(1 into ((42(1 and solving for u>, we find, 



N 



-(J) 



Nrn 

N 



■(J) . 



(48) 



To get an expression for Auj, we consider the energy conservation for the original iV-body model Q. By averaging 
over one complete rotation of the macroparticle, we write: 



N 



(49) 



We then bring into evidence the contributions associated to the massive agglomerate and to the particles of the 
surrounding halo, for both the kinetic and the interaction terms: 



Nm (Placro) , ^ c (p 2 sea ) 



N 



N 2 



2V^7 E 
N . ^ 



sin <?, 



2\/jl E singA =<5(J) 



(50) 



Hereafter we make use of (pmacro) — {Pmacro} 2 , which in turn amounts to assume small oscillations around the 
mean (p ma cro), consistently with l(44|) which neglects such oscillations. The kinetic energy associated to the uniform 
sea can be estimated as follows: 



(P 



2 ) 

sea I 



. A6 / l —dp = ± = - 3w 2 -r 



(51) 



In this estimate we assumed the particles to be distributed uniformly in a rectangular box, disregarding the sinusoidal 
shape of the boundaries. The modulation of the outer frontiers results in higher order corrections which can be 
neglected. 

The interaction term follows directly from l(44() and ((47(1 . Inserting (|51|> in 1(50(1. one finds: 



Au 2 = ^[m{J) 2 ~2A5{J)] 



(52) 
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VI. CLOSED EXPRESSIONS FOR THE AMPLITUDES V± OF THE OUTER BOUNDARIES 

The preceding calculations lead to three equations that allow to estimate N m , uj and Aw (and thus u)+ and 
provided an expression for V± is given. 

Note that V± enter as variables in the Hamiltonian and thus evolve self-consistently. However, one can assess their 
average value using an adiabatic argument £| . The boundaries of the sea are rather sharply drawn by the motion of 
particles, which move following the time-dependent Hamiltonian Hi 5 = p 2 /2 + 2\f~J sm(d — ip), with 1.5 degrees of 
freedom. In a first approximation, the time dependence of J and <p results in a detuning of the wave, as 

ip = djH =-5+ {Niy 1 ' 2 sin^- - ip) = -5 - b/VJ = -J (53) 

3 

where we used (|31|) and the estimates of sec. II VI The resulting velocity agrees with the canonical estimate of reference 
in the low-temperature regime, as it must since ensembles are equivalent for this model. 

Let us neglect the pulsations of J and further variations of (p. Then the Hamiltonian Hi 5 is brought to the 
intcgrable form of the classical pendulum, by a Galileo transformation to the frame moving at velocity v = — (J). In 
this frame, the motion of particles preserves their action, which is directly related to their effective energy 

where p' ± = p± + (J), q' = 6 + (J)t — ip' ; here ip' = ip + (J)t is the wave phase in the new frame. Solving for p, we 
get: 



p± + (J) = ±s/2EZ\^' ^^sing'~±V2£l^l 



smg 



E± 



- (u}± + (J))T\V±\smq' = ±(^-\V±\sinq') (55) 

where we used the assumption that (J) 1 ' 2 /E± is small, i.e. that the particle motion near the sea boundary (which 
has period ~ 2n/ yj2E±) is fast with respect to the characteristic frequency \/2{J) 1 / i of particle oscillations inside 
the wave potential well. From 155fl we finally obtain 

liy-igZ. (so 

Note that the boundaries found here are also approximately level lines of the distribution / given by cq. I|29() . as 
they should be. 

VII. VALIDATING THE THEORY THROUGH DIRECT SIMULATION 

In conclusion, we obtained a system of equations that can be solved to compute the values of V+ and V- and 
provide self-consistent estimates of the three quantities N m , uj + and w_, through a direct calculation of u> and Au>, 
for an initially monokinetic beam and infinitesimal wave {H/N = 0, P/N = 0). The final equations read in the case 
6 = 0: 

Nrn 6-j7(2Ag) 

N l-V7(2Aw) 
Q = -(J) (57) 

\v±\ 



Alu 

V = \V+\ + \V-\ 
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FIG. 6: Rescaled laser intensity I/N vs time. The solid line refers to the reduced dynamics. The free parameters are fixed 
according to the self-consistent derivation here discussed. The dashed line refers to a direct integration of the original N-body 
system 

where (J) = |j4| 2 and b are calculated from l|31|) . The system is explicit after a few algebraic calculations, which show 
that 

K = h-J™^,. (58) 
N 2 V 4 V ; 

It is worth recalling that the derivation of this system assumes that the dynamical variables like V± and J do not 
vary too much, so that these final estimates yield coefficients ui± which do not depend on time. 

In table [I] the solutions of system (|57|l are compared with the values measured in direct numerical experiments. 

The theoretical values are then used to simulate the reduced dynamics (|19f) . Comparisons with numerical results 
based on the original iV-body 10} model are reported in fig. and display a remarkably good agreement. The 
normalized standard deviations of the intensity in the two cases differ by less than five percent. As an additional 
check, we report in fig. |7| the Fourier transform of the signals displayed in fig. A discrepancy smaller than 1 % 
is observed for the peak values. These results provide an a posteriori validation of the reduced model which can be 
effectively employed to characterize the saturated evolution of the FEL field. 





Theory 


Numerics 


b 


0.530 


0.56 


(J) 


0.655 


0.62 


N m /N 


0.472 


0.5 




-0.655 


-0.50 


Auj 


5.408 


5.10 


\v+\ 


0.599 


0.85 


\V-\ 


0.599 


0.8 



TABLE I: Theoretical estimates from system \51\j (first column) vs. numerical measurements (second column) for S = 0. To 
estimate N m we select the i th particle and calculate the spatial distances from the adjacent neighbors. If such quantities remain 
below a cut-off threshold during the subsequent evolution, the particles are said to belong to the inner cluster. 
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Frequency 



FIG. 7: Fourier transform of the signals reported in fig. \Q The thick solid line refers to the reduced model, while the dashed 
one refers to the original N-body system. 

VIII. CONCLUSION 

In this paper we consider a reduced Hamiltonian formulation to reproduce the saturated regime of a single-pass 
Free Electron Laser around perfect tuning. The model consists of only four degrees of freedom. The particles trapped 
in the resonance give birth to a coherent structure localized in space which is formally treated as an individual macro- 
particle. The remaining particles constitute the so-called chaotic sea and are uniformly distributed in phase space 
between two sharp boundaries that evolve according to the Vlasov equation. 

This minimalistic approach was first derived by Tennyson et al. working in the context of plasma physics and recently 
adapted to the FEL in Q . The original formulation assumes three external parameters, namely the mass of the dense 
core, N m , and the average momenta lu± associated to the outer boundaries, which have been so far phenomcnologically 
adjusted to reproduce the results of iV-body simulations. The lack of a self-consistent characterization has significantly 
weakened the potential impact of the reduced Hamiltonian approach so far. 

This long standing problem 0,0 is here addressed and solved. To this end we make use of the analytical estimates 
of the average intensity and bunching parameter at the saturation based on a statistical mechanics treatment of the 
Vlasov equation and of adiabatic theory ["1 11 ll| . Our theoretical predictions for the parameters are compared to direct 
numerical measurements, showing an excellent agreement. Simulations of the reduced dynamics, complemented with 
the strategy here discussed, reproduce remarkably well the oscillatory regime displayed by the original TV-body model. 

In conclusion, we have provided a solid self-consistent formulation for the reduced Hamiltonian, which we hereby 
propose as an alternative theoretical tool to investigate the saturated regime of a FEL. 
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Appendix 

Consider the following expression: 

W = — ^= V sin^ (59) 

i^macro 
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and suppose we want to evaluate (W). The asymptotic intensity oscillates around a given plateau (J), hence we write 
the dominant Fourier term (see fig. 0): 



J = (J) (1 + £,/ sin(wz)) (60) 

with ej = A J I (J), noting A J the amplitude of the oscillations. Moreover the macroparticle displays regular oscilla- 
tions, centered around q = 3n/2. This in turn allows us to write: 



Qi = q + £ Q sin(cjz) (61) 

for the i th particle composing the dense core. Here e q measures the average distance from q. Assuming both ej and 
e q to be small, one can approximate 1)59(1 as follows: 



where use has been made of the fact the cosg = 0. Averaging 1(62(1 over one period z ro t = 2tt/ui one gets: 



N 

(W) = ^= (63) 

N^(jj 



which is identical to equation (|4*5|l . 
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